UTPT-01-02 
hep-ph/0101190 



Renormalon Analysis of Heavy-Light Exclusive B Decays 



Craig N. Burrell and Alexander R. Williamson 

Department of Physics, University of Toronto 
60 St. George Street, Toronto, Ontario, Canada M5S 1A7 



Abstract 

We study two-body exclusive decays of the form B — * D^L (L = ir, p) in 
the heavy-quark limit. We perform a renormalon analysis of such processes 
to determine the order at which nonperturbative factorization-breaking power 
corrections enter the amplitude. We find that a class of leading power cor- 
rections to the color octet matrix element, of ©(Aqcd/ 77 ^)) vanish in the 
limit of a symmetric light meson parton distribution function. We discuss the 
phenomenological significance of this result. 
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I. INTRODUCTION 



The weak decays of B mesons into hadronic final states are important for an under- 
standing of the CKM sector of the standard model, and particularly for the study of CP 
violation. These decays involve a mixture of calculable weak physics, perturbative QCD, and 
nonperturbative QCD. It is the latter component which contains the bulk of the theoretical 
difficulty, and it is primarily contained in the evaluation of low energy matrix elements of 
quark operators. 

A common assumption used to simplify the nonperturbative component of hadronic 
decays is factorization, which specifies that the hadronic matrix element of a four-quark 
operator be factored into two matrix elements of simple currents. For example, 

{it-D+\(&) v _ A {du)v- A \B) - (<K~\(du) v _ A \0)(D + \(cb) v _ A \B). (1.1) 

This prescription (which, following the authors of [0,0, we refer to as 'naive factorization') 
considerably simplifies matters because the factored structures on the right hand side may 
be parameterized in terms of decay constants and form factors. It amounts, however, to 
ignoring corrections which connect the (7r~) to the (BD + ) system. Since these corrections 
are responsible for final-state rescattering and strong interaction phase shifts, leaving them 
out ignores important physics. Also, the left hand side of ( |1 . 1| ) is renormalization scale 
dependent, while the right-hand side is not — a clear indication that relevant physics is 
being lost. 

Recently it was argued that, for certain B decays to heavy-light final states B — > HL, 
the strong interactions which break factorization are hard in the heavy quark limit 
and can therefore be calculated perturbatively. This proposal has been explicitly verified to 
two- loop order ||. This idea allows one to include perturbative corrections missing in (p..l|) 
without introducing any new nonperturbative parameters. It thus provides a remarkably 
attractive means to study rescattering and strong phases in two-body hadronic decays. A 
generalization of this idea has also been proposed for decays to two light mesons HQ], but 
in this article we will restrict our attention to final states with one heavy meson (D, D*) 
and one light meson (7r,p). 

A potential problem with this proposal is that it is valid only in the strict heavy quark 
limit Q. It receives power corrections of the form (AgcD/ 77 ^)™ from a variety of sources: 
examples are hard spectator interactions, non-factorizable soft and collinear gluon exchange, 
and transverse momenta of quarks in the light meson. Unlike power corrections in inclusive 
B decays, there is as yet no systematic way to compute these corrections for exclusive decays. 
By naive power counting one expects such corrections at C^Aqcd/^), but situations are 
known where the naively expected corrections vanish. For instance, in the zero recoil B — ► D 
transition matrix element in Heavy Quark Effective Theory the leading l/m& corrections 
vanish ||. Therefore rather than relying on the naive expectation, one would like to calculate 
the power corrections directly. 



1 Following the authors of [||J2|, we assume that the physical b quark mass is not so large that 
Sudakov form factors modify the power counting. 
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In the absence of a general theory of power corrections, we aim to determine at least the 
order at which power corrections enter. This may be done by carrying out a renormalon 
analysis, which involves calculating a subset of Feynman graphs at each order in perturbation 
theory. It exploits the fact that the perturbative series in quantum field theory is asymptotic, 
and permits one to extract from the large-order behavior of the theory information about the 
scaling behavior of nonperturbative power corrections. In this paper, we use renormalons 
to assess the parametric size of a subset of power corrections to the separation ( |1 . 1| ) for the 
class of B — > HL decays discussed in 0. 

The remainder of this paper is organized as follows: in Section [TI| we review the theoretical 
framework of renormalon analyses. In Section [TTT] we discuss the phenomenological context 
and motivation of our study. In Section [Iy| we describe the renormalon calculation and give 
the result. Section [V] contains our conclusions. 



II. RENORMALONS 

A typical amplitude in QCD perturbation theory may be expressed in the general form 

oo 

R{a s ) = Rn< +1 - (2-1) 

n=0 

Normally one only calculates a few terms in this series. However, one may ask about the 
general behavior of this series at large orders in perturbation theory. It has been argued 0] 
that quantum field theories of phenomenological interest have large order coefficients of the 
form 

R n a n n\ n b (2.2) 

for some constants a and h. Clearly, such a series is factorially divergent. It might appear 
that, as a result, a sum for the series cannot be defined. However, one may define the Borel 
sum R in the following way. Perform a Borel transformation on the series: 

oo oo p 

R{a s ) = £ R n a n s +1 =► B[R}(t) = £ ^t n . (2.3) 

n=0 n=0 n - 

This series in terms of the Borel parameter t is convergent and may be explicitly summed. 
One may then perform an inverse Borel transformation to obtain the Borel sum 

r°° i 

R = / dt e- t/as B\R](t). (2.4) 
Jo 

The original series R and the Borel sum R have the same series expansion. 

In some cases, the transformed series B[R](t) has poles along the positive real axis [|J. 
When these poles are encountered in the inverse Borel transformation ( p.4|) , one is forced to 
deform the integration contour either above or below the real axis. Nothing specifies which 
choice to make, yet the result of the integration depends on the choice. As a result, the 
Borel sum acquires an ambiguity. For a simple pole located at t Q > 0, the ambiguity is 



3 



5R ~ e - to/asM 

(a \ -Wok) /a \ 2m 

where we have defined u = —f3ot > 0. The ambiguity has the form of a nonperturbative 
power correction. For physical quantities, which cannot be ambiguous, there must be present 
power corrections to remove this ambiguity. Note that for // > Aqcd it is the pole nearest 
the origin which gives the leading power correction. This simple sketch illustrates how the 
large-order perturbative behavior of the theory reveals something about the nonperturbative 
sector of the theory ||. 

In general, this approach only permits one to determine the order of the power corrections 
and not their coefficients or analytic form. Furthermore, while a pole at Uq in the Borel plane 
definitely indicates the presence of power corrections ~ (Aqqd / mb) 2u ° , the absence of a pole 
does not necessarily imply the absence of power corrections of that order. The absence of a 
pole is, rather, suggestive that power corrections of the corresponding order are absent [§J. 
The renormalon technique has been applied in a variety of contexts where a general theory 
of power corrections has not been available [[5| p~Q[| . In Section [IV] we study power corrections 
to factorization in this way. 

In practice one cannot sum the entire perturbative series fl2.1| ) to obtain an exact expres- 
sion for the Borel transformed amplitude ( |2.3|) . Instead, one sums a subset of the Feynman 
diagrams at each order, implicitly taking the resulting analytic structure to be characteristic 
of the full result. Typically, all-orders contributions are obtained by inserting into graphs 
'bubble chain' propagators of the kind shown in Figure 1. One may take the formal limit 
Nf — > oo with a s Nf fixed, in which case the set of graphs with a single 'bubble chain' 
insertion are dominant |§. This 'bubble chain' propagator, which we denote by a dashed 
gluon line, has the form |7.JJ 




D,u(k) = 7T2 -9r + -if T,(^N f a s r(H-k 2 /^) + C) n . (2.6) 



n=0 



The factors in the sum arise from the fermion loops depicted in Figure |TJ. These loops have 
been renormalized in an MS-like scheme, and C is a scheme dependent constant. In the MS 
scheme, C = —5/3. We discuss renormalization of amplitudes containing the renormalon 



propagator in more detail in Section [TV]. The Borel transform of this propagator with respect 
to a s Nf is 

B[D^(k)](u) = ( - g , v + E ^(H-k 2 /^ 2 ) + cy 

-(k^k u -k 2 g^). (2.7) 

The limit u — ► of this expression, equivalent to retaining only the first term in the expansion 
depicted in Figure 1, reduces to the usual gluon propagator as expected. For Feynman graphs 
in which the a s dependence arises only from gluon exchange, the Borel transformed graph 
is obtained by replacing the gluon propagator by this renormalon propagator. 
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FIG. 1. The renormalon 'bubble chain' propagator. 

III. FACTORIZATION AND TWO-BODY HADRONIC B DECAYS 

We work in an effective theory where the weak bosons and top quark have been integrated 
out. The relevant part of the effective Hamiltonian, valid below Mw, is 

We// = ^OfeCdOa + CsOs) (3.1) 
where the singlet and octet operators are, respectively, 

Ox = cr M (l - 75 )& d 7 "(l - 7sK (3.2) 
O s = c Tm (1 - 75 )T a 6 d 7 "(l - 75)TV (3.3) 

The running of the Wilson coefficients Cj has been calculated at next-to-leading order |I2| , [13| . 

The computation of any Feynman amplitude in this theory involves evaluating matrix 
elements of the operators 0\$. In general such matrix elements contain both hard and 
soft physics. One would like to disentangle these energy scales; the hard physics could be 
calculated directly and the soft physics parameterized by form factors and decay constants. 
Though it is not obvious that the physics can be disentangled in this way, it has been argued 
that in certain situations it is possible to do so. 

More specifically, consider a decay of the form B — > HL, where if is a heavy meson 
(H = D, D*) and L is a light meson (L = it, p). We require that the topology of the decay 
be such that the light quark in the initial state is transferred to the heavy final state meson. 
In this case, and in the heavy quark limit, it is argued that 'non-factorizable' corrections are 
perturbative and may be calculated. These statements are summarized in the factorization 
equation |l],|2| 

(HL\O l \B) = F B - H (m 2 L )f L [' dxTl{x)® L {x) + ■ ■ ■ (3.4) 

Jo 

where the ellipsis denotes contributions suppressed by powers of AQCD/^b- In this expres- 
sion, the B decay form factor F B ~* H and the light meson decay constant /x are the nonper- 
turbative parameters present in the case of naive factorization (|1 . 1| ) . The 'non-factorizable' 
physics is contained in the convolution of the perturbatively calculable hard-scattering kernel 
T/(x) with the light-cone momentum distribution of the leading Fock state (quark-antiquark) 
of the light meson $^(x). The parameter x is the momentum fraction of one of the quarks 
inside the light meson. 

To leading order in a s , one finds P,[14| 
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I?(x) = l + 0(o2), 



Tl(x) = + O(a a ). 



(3.5) 



Given that the light meson distribution function is normalized to unity, naive factor- 

ization is restored as the leading term in a perturbative expansion. 

There are other decay topologies (penguin, annihilation) for which certain assumptions 
leading to (|3.4j ) are invalid. However, detailed arguments show that these topologies are 
suppressed by powers of Aqcd/^6, and are therefore irrelevant in the heavy quark limit [0. 

In the next section, we present the renormalon analysis of the 'non-factorizable' correc- 
tions. This will involve studying the large-order perturbative properties of (|3.5| ) as a means 
of determining the order at which power corrections enter the factorization equation ( |3.4| ). 

IV. THE CALCULATION 

For a particular heavy final state H (D,D*) and light final state L (tt,p) we must 
calculate the matrix elements 

(0 li8 ) = (H(p')L(q)\0 lj8 \B(p)). (4.1) 



u d u d 




FIG. 2. The factorization-breaking corrections with renormalon propagators. 



The matrix element for the singlet operator is the simplest, so we consider it first. At 
leading order in a s the operator factorizes cleanly into a product of currents. The leading 
factorization-breaking QCD corrections are shown in Figure 0. In order to create a color 
singlet structure from the resulting graph, one would have to consider the light meson to be 
in a Fock state higher than (q'q); this situation, however, is power suppressed by Aqcd/^6 
0. Alternatively, one could retain the leading Fock state and consider the exchange of 
two gluons rather than just one; this situation is suppressed by a s relative to the graphs in 
Figure |2|. In our calculation, then, the singlet matrix element is simply 
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(Ox) = (L(q)\d^(l - 75 )n|0> (F(p')l^(l - 7fe)6|S(p)> 

= ^/L((Jv)-(JA)) (4.2) 

where we define two matrix elements (Jy) = (H(p')\crf ll q ,J, b\B(p)) and (J a) = 
(H(p')\c^q^ 5 b\B(p)). 

The octet matrix element (0%) has a richer structure. It vanishes at leading order, but 
the color structures of the graphs shown in Figure |2| are such that (O s ) receives perturbative 
corrections at all higher orders. Let us denote the individual amplitudes of the diagrams in 
Figure |2| as Ai, and the corresponding Borel transformed amplitudes as -B[^4j]. Defining x to 
be the momentum fraction carried by the up quark of the light meson, x = 1 — x to be the 
momentum fraction of the down quark, and z = m c /mb, the individual Borel transformed 
amplitudes of the diagrams shown in Figure may be written as 



B[A!](u) 
B[A 2 ](u) 
B[A 3 )(u) 
B[A4(u) 



if L C F 


( 


2N C 


\e c m 2 


if L C F 


f " 2 


2N C 


\e c ml 


ihC F 


( f 


2N C 




ifLC F 


{ " 2 


2N C 


\e c m 2 c 



(Fi(x,z,u)(J v ) - Fx(x, -z,u)(J A )) 

{F 2 (x, Z,u)(Jy) - F 2 (x, -Z,u)(Jy)) 

(F 2 (x, l/z,u)(J v ) - F 2 (x, -1/z, u){J A )) 
(Fx(x, 1/z, u) (Jy) - F,(x, -l/z, u)(J A )) 



where 



Fi(x,z,u) 



T(l -2u)T(u- 1) 
T(2-u)(x(l- z 2 )) u 

\-u-2\ 1- 



2m) (3m 



-fi{x,z,u) 



1 



and 



F 2 (x,z,u) 



+ 



r(2 - 2u)T(u - 1) 

r(3-M)(x(i-^ 2 )) u 

2-u ( 2 



x(l — z 1 



2-u + u 2 



+ 



u 



:i+u) 



f 2 (x,z,u) 



2uz 
1 — x(l — z 2 
2uz 



fi(x,z,u) 



1 - 2u \x(l - z 2 ) v 7 l-x{l-z 2 
We have written our result in terms of the hypergeometric functions 

1 \ 



f 2 (x,z,u) }. 



fx(x, z, u) = 2 F X 1 - u, u; 2 - u; 1 



x(l — z 2 



f 2 (x,z,u) 



iFi 1 - u, 1 + u: 2 - u: 1 



x(l — z 2 



(4.3) 



(4.4) 



(4.5) 

(4.6) 
(4.7) 



We have evaluated each of these graphs in d = 4 dimensions, as the divergences which 
would normally be present are regulated by the Borel parameter u in the renormalon prop- 
agator (|2.7|). In the limit u — > 0, however, each of the B[Aj\{u) contain both infrared and 



ultraviolet divergences. It was a central result of [l],FJ] to show that in the heavy quark limit 
the sum 



(4.8) 



i=l 



is infrared finite for amplitudes of the type we are considering. This cancellation may be 
seen explicitly from our amplitudes (fOf): in the vicinity u ~ we have 



B[Ai]{u) 
B[A 2 ](u) 
B[A 3 ](u) 
B[A±]{u) 



if L C f 
2N C 

ihC F 
2N C 

ihC F 
2N C 

ihC* 
2N r 



(Jv-a) 
(Jv-a) 
(Jv-a) 
(Jv-a) 



\ + - (-21ogx(l-z 2 



u 



C + log^-2 

mi , 



- + - 21ogx(l-z 2 ) + C-lo£ 



u 



mi 



- 1 



1 1 



— + - 2 logx(l - 1/z 2 ) + C - log ^ - 1 



u- 



u 



mi 



+ - -21ogx(l-l/^)-C + log^-2 



+ ••• 
+ ... 

+ ■ 
+ ■ 



u 



nit 



(4.9) 



where the ellipses denote terms finite as u — > 0, and we introduce the shorthand (Jv-a) — 
(Jv) — (J a)- The collinear divergences cancel in pairs among the even and odd numbered am- 
plitudes, while the remaining soft divergences cancel in the manner prescribed by Bjorken's 
color transparency argument PJISfl. 

For the sum of the four amplitudes we find 



B [A)(u ~ 0) = - l I^(J v _ A )^ + finite. 



(4.10) 



This remaining divergence is ultraviolet in origin, and may be renormalized in a manner 
consistent with MS-like subtraction schemes. To this end, we follow the prescription of 
|7|jl6[| by defining a renormalized amplitude 



B[A)(u) = B [A](u) + S A (u) 



(4.11) 



where S^(u) contains a divergence which cancels that in S [^4](w) at the origin but is finite 
elsewhere. This regulating function may be written as 



Sa(u) 



E 



9n 



(4.12) 



where the coefficients g n are the expansion coefficients of another function G(e) = J2 n =o9n£ n 
related to the amplitudes computed in d = 4 — 2e dimensions. We refer the reader to the 
relevant literature for an explanation of this method |7|,[i~6f. For the sum of amplitudes A 
we find 



if L C F 2(1 + e )(l + 2e)(3 + 2e)r(4 + 2e) 
{€) 2N C 3r(l-e)r 2 (2 + e)r(3 + e) {V ~ A) 
ifhC F , 127 „ 



(4.13) 



The first term in this expansion cancels the ultraviolet divergence remaining in ( fl.lOD ; 
B[y4.](w) is infrared and ultraviolet finite. 

In this notation, the Borel transform of the factorization equation ( |3.4j ) is 

B[{0 8 )]{u) = I' dxB[A)(u)<S> L {x). (4.14) 
J o 

Comparing to (|3.4|) , we see that i?[y4](w) is proportional to the Borel transform of the hard 
scattering kernel T^ix). 



A. Borel poles and power corrections 



Recall from ( |2.5|) that it is the pole nearest the origin on the positive real axis in the 
Borel plane which indicates the leading power correction. A pole at the origin would indi- 
cate an 0(1) correction, but the renormalization procedure outlined in the previous section 
ensures that no such pole is present. We find that the first pole is located at u — 1/2, 
corresponding to a power correction of C(AQ CD /m fe ). This is not a very surprising result, 
as power corrections of this order are known to be present from a variety of sources 
However, in the vicinity of this pole the Borel transformed amplitude has the form 



B[(0 8 )](tt~l/2)« 



1 



u - 1/2 Jo 



[x — X) 



[Jv-a) + 



(4.15) 



where x — 1 — x, and the ellipsis denotes nonsingular terms. 

Before interpreting this result, we must specify the form of the light meson momentum 
distribution $l(^)- It is customary to write it as [|T/J 



$l(x) = 6x(l - x) 



l + J2^)C 3 n /2 (2x-l) 



n=l 



(4.16) 



where the Gegenbauer polynomials C% 2 (y) are given by 



cT{y) 



1 d n 
n\ dh n 



1 - 2hy + h 2 



-3/2 



(4.17) 



h=0 



The distribution function $l(o;) is a nonperturbative object for which the Gegenbauer 
moments a^(fi) are unspecified. It is known, however, that a^(fi — > oo) = For 
fi ~ rrib ^> Aqcd, one may take the distribution function to have the asymptotic form 
$°(x) = 6x(l — x) up to power corrections of (9(Aqcd/^6)- In this case the wavefunction 
is symmetric under x — * 1 — x, and the integration over x in ( |4.15| ) vanishes, removing the 
pole. 

More generally, note that the Gegenbauer polynomials C^ 2 (2x — 1) in ( [4.16| ) with n 
even are, like $>° L (x), even under x — > 1 — x, while those with n odd are odd under the same 
replacement. The result of the integration ( |4.15| ) may then be written as 



B[(0 8 )](u ~ 1/2) ex — ^(Jv-a) (ofW + « 3 L (^) + of 0*) + " " ") 



(4.18) 
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Only the asymmetric terms of <&l{%) contribute to the residue. Thus the renormalon 
analysis indicates the presence of 0(Aqcd/^&) factorization-breaking power corrections only 
for asymmetric light meson wavef unctions. 

For certain final states of interest (L = n, p), SU(2) symmetry ensures that the wave- 



function is symmetric |I8|.|T9f. In these cases the pole at u = 1/2 disappears, and the 



renormalon analysis gives no indication of factorization-breaking 0(Aqcd/^&) corrections, 
suggesting that such power corrections are likely absent. This result implies that power 
corrections due to the 'non-factorizable' vertex diagrams in Figure |2| should be suppressed 
in B — ► D^ + ix~ and B — > D^ + p~ decays relative to decays into a light meson with an 
asymmetric distribution function, such as B — > D^ + K~ . 

In assessing this result, one should keep in mind that there are other sources of power 
corrections to these decays which our analysis does not address. For example, we have not 
considered soft interactions with the spectator quark, annihilation diagrams, or interactions 
with sub-leading Fock states involving additional soft partons |T],0. Our conclusions about 
the vanishing C^Aqcd/^) corrections in the symmetric limit apply only to the class of 
diagrams we have considered. 

The next renormalon pole is located at u — 1, corresponding to an 0(Aq CD /r??|) cor- 
rection, and here there is a nonzero residue even in the case of a symmetric wavefunction 
$l{x). Therefore one should expect power corrections of this order for all decay modes. 



V. CONCLUSIONS 

In this paper we have carried out a renormalon analysis of factorization-breaking effects 
in B meson decays to certain hadronic heavy-light final states. The renormalon approach, 
which probes the theory at high orders in perturbation theory, allows us to learn about 
nonperturbative power corrections. In the low energy effective theory governing the B 
decays there are two operators, a color singlet and a color octet. The factorization-breaking 
corrections to the singlet operator are, however, suppressed by powers of ot s or AqcD/^fe 
relative to the octet operator; for this reason we focus our analysis on the octet. 

We find that the renormalon analysis of 'non-factorizable' corrections to the octet matrix 
element indicates the presence of power corrections of C(Aqcd/^&)- We also find that these 
leading power corrections are sensitive only to asymmetries in the light meson light-cone 
parton distribution function. Thus for a symmetric distribution function the leading renor- 
malon pole vanishes, suggesting the absence of (9(Aqcd/^&) 'non-factorizable' corrections 
in such cases. The next power corrections are present at 0(Aq CD /7??|). 

As it is natural to expect the light meson wave function to be symmetric for certain 
final states (L = ir,p), the potential vanishing of the 0(AQCT>/mb) power corrections in the 
symmetric limit is an interesting result that warrants a few additional comments. First, it 
should be noted that a particular decay B — > HL of the type we have considered in this 
paper always involves a combination of the singlet and octet operator matrix elements. Our 
result, however, applies only to the octet matrix element — we are unable to conclude from 
our analysis whether the leading power corrections to the singlet matrix element vanish in 
an analogous way. 
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Second, the power corrections we can probe through the renormalon analysis are phys- 
ically due to soft and collinear 'non-factorizable' gluon exchange, and our comments about 
suppression of the leading corrections should be understood as referring to these effects only. 
There are, however, other 0(Aq CD /m fe ) corrections originating from diagrams which were 
not considered in our analysis. A more detailed examination of these effects would have to 
be undertaken to accurately assess the size of ©(AqcdA 7 ^) corrections to a physical decay 
amplitude. 
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